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Abstract 

We perform a phenomenological analysis of the decays B — > D**n, where 
D** is a P-wave excited meson with total angular momentum j = ^ or | for 
the light cloud, recently measured by the Belle Collaboration in the modes 
B° -> D** + n- (Class I) and B~ -> D** -k~ (Class III). Making the reason- 
able assumption of naive factorization, that we test in B — » D(D*)tt decays, 
Class I decays allow to extract the Isgur-Wise form factors Tx/ 2 (w), t 3 / 2 (^) at 
w = w max (q 2 = 0). We obtain T 1/2 («Wr) < 0.20, T 3 / 2 (w max ) = 0.31 ± 0.12. 
We discuss the question of the w dependence of these IW functions. We find 
agreement with the Bakamjian-Thomas quark model of form factors and, ex- 
trapolating at w = 1, with Bjorken and Uraltsev sum rules. We discuss also 
Class III decays, where the D**° (j = |) emission diagram contributes. We 
extract the corresponding /d 1/2 decay constant, that is in agreement with 
theoretical estimates at finite mass. Finally, we must warn that 1/ttiq correc- 
tions could be large and upset the results of the present stage of this analysis. 
On the other hand, we confront present data on the semileptonic rate of B 
mesons to excited states with theoretical expectations. 
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1 Introduction. 



Our main purpose in this paper is to extract information on the lowest (n = 
0) heavy quark Isgur-Wise functions Ti/ 2 (1) and r 3 / 2 (l), that correspond to the 
transitions ~ — > | + or | + , quantities that are of importance in heavy quark physics, 
for instance in small velocity sum rules (SR). We will first make a naive estimation 
of these quantities from B° — > D** + tt~ decays assuming factorization and the heavy 
quark limit. We will moreover discuss the more involved decays B~ — ► D**°tt~ in 
which there is D**° emission, in particular the sign of the interference between the 
7r~ and D**° emission diagrams. 

There are four P-wave D** mesons corresponding to the coupling of the total 
quark spin S = 0, 1 and the orbital momentum £ — 1. In the language of the 
heavy quark limit, where the total angular momentum j of the light quark is a good 
quantum number (j = ~ or |), these states can be denoted by Dj where J is the total 
angular momentum of the state. There are then four possibilities (j,J) = (§;0), 
(1,1), (|,l) or (|,2). 

According to the states of charge, the B — > D**n decays are of three classes, 
following the classification of B. Stech and collaborators pQ, [2] : 

(i) B° — > D** + 7i~ (Class I) where only the 7r emission diagram contributes ; 

(ii) 5° — > D**°7i° (Class II) where only the D**° diagram contributes ; 

(iii) B~ — > D**°7T~ (Class III), where both diagrams of 7r emission and D**° 
emission contribute. 

It is worth to recall that isospin symmetry relates the amplitudes of the three 
classes in this particular case, namely [S] : 

A (b~ -> D**°ti-) = A (S° -> D** + n-) - V2A (B° -> D**°n°) . (1) 

The Belle Collaboration has obtained in the past very interesting results on 
Class III decays B~ — ► D**°tt~ 4J. Four states where indeed observed, two narrow 
states corresponding to the j = | states, that decay into D(D*) in the D-w&ve, and 
two very wide states that decay into D(D*) in the S'-wave. 

The interesting news is that recently, at the Beijing ICHEP 04 Conference, the 
Belle Collaboration has presented results on Class I decays B° — > D** + n~ [5J. In- 
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terestingly, the narrow states have been observed, while only limits on the decays 
into wide states have been obtained, indicating a much smaller BR. The Belle data 
have recently drawed the attention of the theory [Hj. 

An enormous theoretical effort has been dedicated in the last five years to the 
understanding of non-leptonic two-body B decays in the cases of the emission of 
a light meson like it or p in the so-called QCD Factorization approach [Jj, in the 
Perturbative QCD Factorization approach [8 or within the Soft Collinear Effective 
Theory |9 a . These methods have been applied to two-body decays into ground state 
mesons. In the present paper we are dealing with decays into excited D** mesons, 
with both 7r or D** emission diagrams. For Class I decays we have only the 7r 
emission diagram, and in this case we could in principle use the QCD methods of 
these papers. However, in Class III decays there is the diagram of D** emission, 
a meson composed of heavy-light quarks, for which there are no rigorous results. 
Moreover, we are dealing with the first measurements of these decays, that hopefully 
will be refined in the future. For these reasons, as a preliminary study, we will stick 
to the naive factorization approach pQ, j2] in order to investigate if there is a sensible 
description of the decays B — > D**tt within this simple phenomenological approach. 

Class I decays and the 7r emission diagram of Class III are related to the B — ► 
D** form factors that, in the heavy quark limit, reduce to two Isgur-Wise (IW) 
functions 71/2(10) and t$/2{w) |10j . These form factors are of a significant theoretical 
importance, since they are related, at zero recoil w = 1, to the slope of the elastic 
IW function through Bjorken [TUj . [TT] and Uraltsev SR |12j . 

In a recent paper we have tried to use the Belle data on Class III decays to extract 
T~i/2(w max ) and T 3 /2(w max ), where w(q 2 = 0) = w max [TT?| . This calculation relied 
on a strong hypothesis, namely that the diagram of D**° emission should be small. 
We got some results on Ti/ 2 (w max ), T 3 / 2 (w max ) that were extrapolated to 71/2(1), 
r 3/2(l) assuming the w-dependence of the form factors given by the Bakamjian- 
Thomas (BT) class of relativistic quark models that yield covariant form factors 
exhibiting heavy quark symmetry IT II. We obtained Ti/ 2 (1) ~ r 3 / 2 (l), at odds with 
the expectations of Uraltsev SR. 

However, to neglect the D**° diagram is a rough approximation that could be 
unfounded ^H]- It is well-known that in some cases the color-suppressed diagrams 
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like the D**° emission one are often not as suppressed as one could expect on naive 
grounds. Therefore, our determination of Ti/ 2 (w max ), T 3 / 2 (w max ) has to be recon- 
sidered using only Class I decays, now measured, and where only the 7r emission 
diagram contributes. Hence the interest of the new results on Class I decays on 
which we will first concentrate. Below we will come back to the interpretation of 
the results on Class III decays, taking into account D**° emission. In what follows 
there is some unavoidable overlap with our Appendix B of ref. ^3] from which some 
points are worth to be recalled. 

The paper is organized as follows. In Section 2 we extract Ti/2(w max ), T^/2{w max ) 
from Class I decays 5° — > D** + n~ using factorization, and discuss the question of 
the ^-dependence and the extrapolation of ri/ 2 (w), t 3 / 2 (u>) at w — 1, the comparison 
with Bjorken and Uraltsev SR and with BT quark models. In Section 3 we combine 
the results of Section 2 with the measured rates of Class III decays B~ — > D**°tt~ 
in order to extract the needed value for the decay constant /d 1/2 , non- vanishing 
in the heavy quark limit, and we compare with theoretical predictions. We treat 
with special care the question of the interference between the tc~ and D**° emission 
diagrams. In Section 4 we make predictions for Class II decays B° — > D** 7t°. In 
Section 5 we discuss the implications for the semileptonic decays B° D** + £~V and 
compare with the scarce existing data, and in Section 6 we conclude. In Appendix 
A we reproduce the Belle data for Class I and Class III decays to make clear in the 
text how we extract the rates B — > D**n of the different modes and how we treat the 
experimental errors. In Appendix B we use the factorization model to compare with 
the data on B — > D(D*)tt, where all modes have been measured, and we extract the 
effective coefficients a\ and a 2 that enter in the factorization model. In Appendix 
C we discuss the corrections to factorization and finally Appendix D is devoted to 
the question of the l/mg corrections to the heavy quark limit. 



decays. 

Let us now consider the Class I decays measured by the Belle Collaboration, 
where for the wide states the masses of the results of Class III decays are assumed 



2 Extraction 




from Class I 
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(Table 5 of Appendix A). 

Assuming that these states decay essentially into two-body modes, i.e. 
B (Df 2 -> (D + D*)ir\ B (l> 3/2 -> L>*tt), £ (Dq /2 -> dA 5 (d\ /2 -> £>*7r), the 
following branching ratios are given by a Clebsch-Gordan coefficient 

fi (£,3/2 + ^ ^O^) = ^ ^1/2 + ^ = B ( D l/2 + ^ ^O^) = 2 (2) 

To estimate B (Df 2 + -> £>°7r+) and B (d 3 / 2 + -> D*°7r + ), we use the spin 
counting of the non-relativistic quark model. In the limit of heavy quark symme- 
try, i.e. assuming the pairs (D,D*), (D^ 2 , D^ 2 ) and (D^ 2 , D^ 2 ) to be degenerate, 
simple angular momentum calculations give 

r (pl /2 ) = r (d 3/2 ) r (dI /2 ) = r (z^ 2 ) (3) 

r (£>f _> ZTtt) = I r (D 2 3/2 - Dtt) . (4) 

This last relation gives the needed spin counting coefficient. 

It is easy to obtain this factor by realizing that to have the D wave one needs 
(1 denoting the quark emitting a pion) the operator (taking Oz along the pion 
momentum) to emit a pion in the D wave reads (erf k w ) exp(i^iA; 7r ) — > ik^afzi. We 
have then, for the non-vanishing amplitudes 

M (£>2 /2 -> £>7r) = < 1 0, 1 0|2 > < 1 erf 1 1 > < 0|^f|l > 
M / jD 3/2(±i) _^ £,^±1)^ =< i 0, 1±1|2±1 > < l±l|af|l±l > < 0|X|10 > (5) 
that gives 

m (pT -+ Dir) = y| < 0|^r|l > 

M (dI /2{±1) -> ZT^tt) = < 0|^|1 > (6) 

and hence (jlj). 

We now take into account the actual masses. Since both D^ 2 — > Dir and -D 3 / 2 ~~ * 
D*ir proceed through the .D-wave, we will have, in an obvious notation, in the isospin 
symmetry limit, 



r 


( D r + 


-> 1^+) 


r 


(pT - 


-> D + 7T-) 


2 


p5 


r( 




-> Z> 7T+) 


r( 


^3/2 _ 


-> L>*+A 


! 3 





= 2.5 . (7) 
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This estimation is in agreement with the present world averages fT6J 

1.9 ± 1.1 ±0.3 

2.3 ±0.6 . 
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( D r + 


-> D°Tl + ) 


r 




-> D*°TT + ) 



YD 



3/2 



D~ 



71 



r(D 



3/2 



7T 



(8) 



Therefore, we obtain the branching ratios 



B D 



)3/2 



71 



0.48 



,3/2 



71 



From these BR, adding the errors in quadrature, we find roughly 



3/2 



3/2 



0.19 . 

— ► D°7l 

1*0 



(9) 



B (b° -> Dl' 2 + tt-) = (6.4 ± 0.8) x 10" 4 (from D 

B (B° -> L>2 /2 + tt-) = (12.9 ± 3.3) x 10- 4 (from D 

We realize that the B (5° — > D^ 2 +7l ^j differs if one obtains it from D 



Zr°vr + (10) 



3/2 



D°7i + or from D 



3/2 



D*°7r + , and the values are consistent only within 2a. Using 



PJ) for the other modes, and taking into account the large uncertainty from both 
results (jTUJl one finds 











(b°- 






(b°- 


-> D 1/2 




(b°- 





71 

7T~ 
7T~ 



= (10.9 ±5.3) x 10~ 4 
= (5.5 ± 1.4) x 10~ 4 

< 1.8 x 10~ 4 

< 1.0 x 10~ 4 . 



Assuming factorization of tt~ emission, as it is reasonable within the BBNS QCD 
factorization scheme in the heavy quark limit [7^ and assuming that the states 1 + 
are unmixed, we find for the decay rates 

r = ^ 2p 

16n nr, 



B 



12) 



where 

m(b° 



M (B° 



3/2 



M [B° 



Di 



1/2 



M (~B° 



D 



1/2 



71 
7T~ 



2m D m B (m B + m D ) 2 (w 2 , - l) a\fl 
2m D m B (m B - m D f (w + 1)' 
4m D m B (m B - m D ) 2 (w% - l) a\fl 



T3/2 (w ; 



2 

7T 

2 „2 r2 



o 2 " 1) a?/. 2 



T"3/2l^0j 



Am D m B (m B + m D ) (w - 1) t 1/2 (w ) 
6 



(13) 



with 



m 2 B + m 2 D 



= -r ( 14 ) 

the subindex denoting the value of w for q 2 = m 2 = and mjj the mass of the 
corresponding Dj state. The short distance QCD factor a x is close to 1 (Appendix 
B). 

It is interesting to notice that the rates p2j) . (fT3*|) are given, assuming the D** 
for a given j = \ or § to be degenerate, by the expressions 



G 



16tt 



D 



3/2 



7T 



r [B 



F it/ t/* I 2 3 2 ^2 (1 



n 3/2 + - 

l 5 (l + r) 7 



16r 3 



T3/2 



1 + r 2 

2r 



(15) 



r(# 



n l/2 + - 



r(5 



^ W V*\ 2 m* a 2 f2 (l-r) 5 (l + ^ 



2r 



Tl/2 



where r 
equalities V [B 



3/2 

—2- and r 



1/2 



1 + r z 
2r 



(16) 



respectively in the first and second relations. The 



n 3/2 
£>2 



-- r (5 -> dI /2 tt), t(b^d\ /2 <k) = r(s 

follow from heavy quark symmetry since the .B meson is spinless, and there is a 
single helicity amplitude for each decay, the emission of a longitudinally polarized 

Using the central values for the masses, but taking into account the errors in 
(fTTj) and \Vcb\ = 0.040 ± 0.002 (Appendix B), we find from the different modes, 



B 
B 
B 
B 



n 3/2 

Dl'\ 
A) 7T 



73/2(1.31) 
7-3/2(1-32) 

n/ 2 (1.37) 

71/2(1.32) 



= 0.32 ±0.10 
= 0.23 ±0.04 

< 0.20 

< 0.16 . 



(17) 



Within la there is consistency between the different determinations of |t3/2(wo)| 
and |ti/ 2 (wo)|; but errors increase considering both determinations. Since besides 
the statistical errors, there are systematic errors (from experiment and theory), we 
consider safer to take the union of the domains (|17|) rather than their intersection. 
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We conclude safely that we will have the numbers 

r 3/2 (1.31)| = 0.31 ±0.12 
7i/2(1.37) < 0.20 . 



2.1 Extrapolation at w = 1 and comparison with Bjorken 
and Uraltsev sum rules. 

Our results have been obtained at w max . It is of interest to know the values of 
71/2(10)) T3/2 (u>) at w = 1. There are a number of values for the slopes of these IW 
functions in the literature. 

We should keep in mind that we have two rather loose constraints on 71/2(10), 
73/2(10), namely the values at w max (fTHj) and the qualitative idea that the n = IW 
functions should give a main contribution to Bjorken and Uraltsev sum rules. 

Let us consider, as an illustration, the parametrization obtained within BT quark 
models (last reference [T3] ) 

7-1/2(10) =Ti/ 2 (l) 
73/2(l0) = 7-3/2(1) ( 

we obtain at zero recoil 



2 \ 2f7 3/2 
W+l ) 



a( /2 = 0.83 



cr 3/2 = 1.5 



(19) 



73/2(1) 

n/ 2 (i) 



= 0.46 ±0.18 
< 0.26 



(20) 



to be compared with the values in the BT model 

BT 



r 3/2(l) 

n /2 (i) 



BT 



0.54 



0.22 . 



(21) 



We find agreement for |t 3 / 2 (1)| within errors, and |ti/ 2 (1)| could still be roughly 
consistent with the BT model. 

Bjorken JU], JT] and Uraltsev sum rules write, respectively 



P 



(n) 



i + EK/2( 1 )| + 2 E T3/2 



£h$(i 



n 
|2 



Ek ( /1(i)| 



1 

I 



(22) 



It is understood that these SR are truncated at some n that corresponds to a 
scale A ~ 1 GeV and it is then natural to assume that the ground state dominates. 
Keeping thus the n = states, with which we are dealing here, we get contributions 
to Bjorken and Uraltsev SR that lie in the following ranges : 



0.40 < - + 

4 



2 



r 1/2 (l) +2r 3/2 (l) 



2 



< 1.11 



0.03 < 



2 



2 



73/2(1) - r V2 (l) < 0.41 . (23) 



Therefore, the n = states could give an important contribution to the SR and, 
considering this piece as dominant, low values for p 2 are not excluded nor the value 
| for the r.h.s. of Uraltsev SR. 

There are other theoretical estimates of the IW functions ri/ 2 (w), t 3 / 2 (w), mainly 
within the QCD Sum Rules approach [21], [2H], |2E], |2Zj- The pioneering calcula- 
tions of Tyz(w) and r 3 / 2 (w) [2H| show indeed that at large w the slope of Ti/ 2 (w) 
is much smaller than the one of 73/2(117), as in the BT model, while the values at 
w = 1 were found roughly equal for both IW functions, 7i/ 2 (l) = t 3 / 2 (1) = 0.24. On 
the other hand, next-to-leading calculations of the function rx/ 2 (io) have later been 
performed [2E], giving 7i/ 2 (l) = 0.35±0.10 and a slope of the order crf, 2 = 0.5, in our 
notation ()19j). In view of the importance of the corrections, it would be interesting 
to have the corresponding calculation for r 3 / 2 . This latter value for 7~i/ 2 (l) is larger 
than the value obtained in the present paper. On the other hand, in ref. j22j there is 
a calculation of both 7i/ 2 (w) and r 3 / 2 (w) (((w) = 2ti/ 2 (ii>), t(w) = V3t 3 / 2 (w) |28j) 
and gives r 3/2 (l) = 0.43±0.08, r 1/2 (l) = 0.13±0.04 and the slopes oJ /2 = 0.90±0.05, 
a i/2 = 0-^0 ^ 0.05. These results imply a sizeable contribution to Uraltsev SR, are 
in agreement with the determinations of the present paper, and are qualitatively 
consistent with the BT model results (|2T|) and with ()19)) . t 3 / 2 (u>) being steeper than 
r i/2(w). One should notice that a different interpolating field for 7i/ 2 is used by j^Zj 
from the one in [23], [2E], and that radiative corrections are absent. Recently, a lat- 
tice determination has obtained the values Ti/ 2 (1) = 0.38(5) and t 3 / 2 (1) = 0.53(8), 
with unknown systematic errors [211! • These values imply a sizeable contribution to 
Uraltsev SR. Compared with the BT determination ()21|) . r 3 / 2 (l) is in fair agreement, 
while 7i/ 2 (l) is larger, and in agreement with the QCDSR determination j2E]- A 
fortiori, this latter value is much larger than the QCDSR result [27], and also than 
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the present phenomenological limit (fTHj) obtained in the present paper, with the 
extrapolation from w max assumed here. We summarize the situation in Table 1. 



Theoretical 
method 


n /2 (i) 


*l/2 


rs/a(l) 


2 

a 3/2 


QCDSR 
(NLO) [H] 


0.35 ±0.10 


0.5 






QCDSR HZ! 


0.13 ±0.04 


0.50 ±0.05 


0.43 ±0.08 


0.90 ±0.05 


BT Quark 
Model [Hj 


0.22 


0.83 


0.54 


1.5 


Lattice 


0.38(5) 




0.53(8) 




Present 
paper 


< 0.26 


0.83 
(input) 


0.46 ±0.18 


1.5 
(input) 



Table 1. The values at zero recoil Tj(l) and slopes crj (j = |, |) in the different 
theoretical approaches, compared with the phenomenological determination at w max 
of present paper, extrapolated at w — 1 with the slopes of the BT quark model (fTU|) . 



3 Interference with D**° emission in Class III 
decays. 

Let us now consider Class III decays measured by the Belle Collaboration, that 
we summarize in Table 6 of Appendix A. 

We have here the same BR as for the modes of Class I, namely 

B (dT - ZT+7T-) = B (Dl' 2 - D+Tr-) = B (D\ 12 - ZT+7T-) = I 
B (Dl /2 -> = 0.48 5 (Ds 72 -> D*+7T~) = 0.19 . (24) 



We find, adding the errors in quadrature 

5 (B- D2 /2 V) = (7-1 ± 1-6) x 10~ 4 (from D 3 2 /2 -> D+tT 

5 (b~ -> D 2 3/2 V) = (9-5 ± 2.5) x 10~ 4 (from Z^ 2 -> ZT+Tr) .(25) 



We realize that the values obtained for B(B — > D^ 2 °7r ) from -D^ 2 ° — > -D 

^3/2 o 



7T 



or 7r agree within lcr. Using for the other modes, and taking 

10 



into account the uncertainty from both results (|25j) one finds 



B 




D'f- V) 


= (8.7 ±3.2) x KT 4 


B 


(B-- 


-> D?' 2 V) 


= (10.2 ±2.3) x 1(T 4 


B 


(B-- 


n l/2 -A 
A) j 


= (9.1 ±2.9) x 10 -4 


B 




n l/2 _> 

->• D l * j 


= (7.5 ± 1.7) x 1(T 4 . 



(26) 

Comparing these BR with the corresponding Class I (|11|) we see that there is a 
large difference for the j = \ states, while there is consistency for the j ' — | states. 
This is interesting and seems to indicate that the D**° emission diagram could be 
very important [Ej. This is likely, because the decay constants of j = | states do 
not vanish, while those of j = | states vanish in the heavy quark limit jTHJ [19J |3 j 

fi/2 ± /a/2 = . (27) 

As demonstrated in [T%] . the equality / 3 / 2 = follows intuitively from the fact 
that the multiplet j = ~ contains two states with J = 1, 2 and there is no current 
coupling the vacuum to J = 2. On the other hand, fx/2 ^ follows because the 
D**(J = 0) is a system of widely unequal masses, and vector current conservation 
does not hold. We assume, following [T3], that the decays B~ — > D^ 2 °7r~ and 
B~ — > D 1 ^ 2 °7T~ have a sizeable contribution from D 1 / 2 ( J = 0, 1) emission via, 
respectively, the vector and axial current. 

We now consider both diagrams for Class III decays and we will take care of the 
delicate question of the relative sign between the 7r emission and the D**° emission 
diagrams. 

The decays B — > D^ 2 tt and B — > D^-k are respectively S'-wave parity con- 
serving and P-wave parity violating. Let us define in an homogeneous way the 
needed matrix elements (q = p — p') JB|- For 71 emission we need the current matrix 
elements 

< n(q)\A^\0 > = f n q^ 
< Dl 12 \p')\A^\B-(p) > = v^^i 2(v' - «)„ rx/ 2 (w) 
< D\ 12 \p',e)\V^\B~{p) > = 2 \(w - - (e* ■ v)^] r 1/2 (w) (28) 



11 



while for D** emission [T%| . [T7j 

<Dl /20 ( P ')\v,\o> = f Dl/2 P ; 



< D\ 12 V,e)|A M |0 > = -/ Dl/2 m^ /2 

< tt-^I^IB-Cp) > = [p, + % - ff( P ' 2 ) + ^^V,fo B {v 2 ) ■ 

(29) 

The minus sign for the definition of the D 1 ^ 2 decay constant comes from the 
Clebsch-Gordan convention of Isgur and Wise (coupling the orbital angular momen- 
tum £ = 1 with the light quark spin s q = 5 to give j — |) that yields the definitions 
From (|2*9*j) . as predicted by heavy quark symmetry, one obtains, 



< Dl' 2 °(p')\A 3 \B-(p) > = -< D\ 12 °(p',e)\V°\B-(p) > = ^bW^ 2v' 3 r 1/2 (w) 

(30) 

and corresponds to the convention 

S c 3 \dI + > = -\Dl > . (31) 

Likewise, one must have 

< D 1 ' 2 V)l^ 3 |0 > = - < D\' 2 °(p',e)\A°\0 > = f Dl/2 p' 3 . (32) 

This is the convention that we have used in ^H] , (there is a misprint in formula 
(14) of [IH], corrected in |T7|). 

We find for the rates (only one helicity amplitude contributes to the B~ — > 
D^ 2 °7r~ transition) : 

r(s -^l> tt \ = — \v cb v ud \ — 

v y lo7r m B 

a ly /m D m B 2 (m B + m D ) (w - l)f v n/ 2 (w) + a 2 m 2 B f Dy2 f^ B (m^) (33) 



Y(B-^D\' 2 \-) = ^\V cb V: d \ 2 P p2 



a ly /m D m B (m B - m D ) / 7r r 1/2 (w) + a 2 m B m D f Dy2 fl B (mj^ 2 . (34) 

We will use in these expressions the color-allowed and color-suppressed factors 
respectively of the order a\ = 1, a 2 = 0.3 (B.8). The powers of p indicate that the 
decays B — ► D\j 2 ix and B — > .D^V occur respectively in the S and P waves. 
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The relative sign between both terms in ()33|) and ()34|) is crucial. Let us give an 
argument that shows that the interference is constructive. Instead of considering 
the n, let us consider the pseudoscalar D meson, composed of heavy-light quarks. 
Our assumption is that the form factors and decay constants between ground state 
mesons do not change sign when going from heavy mesons made of heavy-light 
quarks to light mesons made of equal mass quarks. This is a very sensible continuity 
hypothesis that is satisfied in the quark model, since there are no nodes in these 
ground state wave functions, and the extrapolation in reduced mass is smooth. 
On the other hand, this smooth continuity in mass is commonly used in lattice 
calculations, and it is also observed, considering for example the decay constants 
f Dg or f K and varying the c or s quark masses. 

In ^H] , \H\ we did demonstrate (using duality in B° — B° mixing and also within 
the OPE) the following sum rules in the heavy quark limit of QCD for heavy-light 
form factors and decay constants, valid for all values of w : 

E f in) Z in) H= 2 E /$ A% (w) = /«» (35) 

n n 

where n denotes a radial quantum number, f^ = f is the ground state decay 
constant and ^°\w) = £,(w) the elastic Isgur-Wise function. The decay constants 
and fy\ scale like in particular, we have demonstrated that the rigorous 

SR (|33j) are satisfied within relativistic BT quark models and in the non-relativistic 
quark model. Within BT quark models, we have shown that a main contribution to 
the SR ()35|) comes from the n = states, that has the same sign as the whole sum 
[TTj and the same is true in the non-relativistic quark model : 

Sign[/ £(w)] = Sign [/ 1/2 r l/2 (w)\ (36) 

where we have used the notations Ty\(w) = Ti/ 2 {w), f[°l = f\/ 2 - Multiplying 
the equalities (|3l)j) by £(w)ti/ 2 (w), we have Sign{/ri/ 2 (w)[^(w)] 2 } = Sign{fi/ 2 £(w) 
[ti/ 2 (w)] 2 } and therefore 

Sign [/ r 1/2 (w)\ = Sign [/ 1/2 £(w)\ . (37) 

Heavy quark scaling implies for B —* D form factors : 



VAmDmB ;fV) = VAmDmB /n ? = £M (38) 



m D + m B m D + m B 1 



(m n +m B ) 2 
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and therefore Sign[f BB (q 2 )} = Sign[/^ B (q 2 )] = Sign[^(w)]. Our continuum assump- 
tion linking heavy-light mesons to light mesons implies then, within a definite phase 
convention : 



Sign [fl B {q 2 )] = Sign [f£ B (q 2 ) 
Sign [f D ] = Sign [f v ] . 



Sign [H(w)} 



(39) 



From (J37)) and ((HSj) we get 



Sign f n T 1/2 [ 



w 



Sign f 1/2 fl B {q 2 ) = Sign f 1/2 f^(q 2 ) 



(40) 



Therefore, a relative constructive sign between the two contributions in ()33|) and 
dSU) follows from (1301 . 

We need some input on the form factors fo B {q 2 ) and f+ B {q 2 )- We could use the 
simple theoretically motivated pole-dipole parametrization for fo B (q 2 ), f+ B {q 2 ) of 
the Large Energy Effective Theory (LEET) [201 : 



(41) 



However, there is an empirical parametrization, inspired by (|4*Tj) . that fits the 
lattice data on these form factors, proposed by Becirevic and Kaidalov j2I] : 

c B (1 - ob) 



/?V) 



c s (1 - a B ) 



Ha B(0+) 



(42) 



A fit to the lattice data [22] yields two sets of values for these parameters. We choose 
one of them, the other one yielding very comparable results : 

c B = 0.51(8)(1) 
a B = 0.45(17)«; 06 



Pb = 1.20(13)^; 
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+0.15 



00 



that corresponds to 



/r(0) = / ^(0) = 0.28(6)(5) 



(43) 



(44) 
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Concerning the QCD coefficient a 2 , in Appendix B we have made an analysis of the 
well-measured decays B — > D(D*)tt in all its charged modes. Since the perturbative 
estimation of a 2 = 0.2 appears to give too small Class II branching ratios, we 
have to consider, following [2], [3] non-perturbative contributions to a± and a 2 (see 
the discussion in Appendices B and C), that yield the values 

ai = 1 a 2 = 0.3 . (45) 

Moreover, we adopt, like in Appendix B, the parametrization for the form factors 

fl B (q 2 ), fS B (l 2 ) g iven by @2MEJ- 

Once we know the sign of the interference between the two terms in (}3*3*j) and 
(I34J) . we proceed as follows. We extract the decay constant /d 1/2 from these formulas 
comparing to the Class I ones 

1/2 o__\ G\ |T/ y 1 2 V r - -i 2 



T[B-^Dr* =Tf- V cb V udl - 
' 1d7t m 



B 



aiy/m D m B 2 (m B + m D ) (w - 1)Uti/ 2 (w) 

(46) 

T -> L>i /2 °7T _ ) = — - iKb^l 2 2" \aiy/m D m B (m B - m D ) f w t 1/2 (w)] . 

v 'An m B m D 1 J 

(47) 

Adding the theoretical errors in quadrature and using the QCD coefficients (|43j). we 
find 

f Di/2 = (206 ± 120) MeV (from B -> Dj /2 7r) 

/ Di/a = (196 ± 93) MeV (from 5 -> D\ /2 tj) . (48) 

Both determinations are roughly consistent. In view of the large systematic uncer- 
tainties, we proceed as in (|18)1. taking the union of both domains rather than the 
intersection. We thus keep the safe range 

/01/2 = (206 ± 120) MeV . (49) 



3.1 Comparison with theoretical estimates of fD 1/2 - 

The value ()49|) is in reasonable agreement with the calculation of QCDSR [24J 
that gives, for decay constants of D mesons with 0~ and + quantum numbers, 
including 1/rriQ and a s corrections, the following numbers 

f D = (195 ± 20) MeV / D{0 +) = (170 ± 20) MeV . (50) 
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There are also calculations within QCDSR in the heavy quark limit, without 
including a s corrections (SHj, that give a larger value for fn(o+), consistent within 
la with flU 

f D = (0.21 ± 0.03) GeV 3/2 
y/m5^) /d(o+) = (0.46 ± 0.06) GeV 3/2 . (51) 

Another estimation using QCDSR in the heavy quark limit [SD] gives a larger value, 

fD(o+) = (0.570 ± 0.08) GeV 3/2 (52) 

and correcting for the Btc continuum (30] one gets the results 

f D = 0.35 GeV 3/2 
^rrr^T) f D(0+) = (0.36 ± 0.10) GeV 3/2 (53) 

that are consistent with (|49|) . 

Within the Bakamjian-Thomas class of relativistic quark models [H], the decay 
constants of heavy-light mesons in the heavy quark limit have been computed jT7] . 
One finds heavy quark scaling ^/m^f^, = Const, for the D(D*) states and for the 
doublet of j p = | + states. Within the specific spectroscopy model of Godfrey and 
Isgur pH], one finds the following values for the lowest n = states : 

v /^ f D = (0.670 ± 0.020) GeV 3/2 

fo 1/2 = (0.640 ± 0.020) GeV 3/2 . (54) 

These values for | and | + doublets are close. From the masses m^D*) and 
rriD 1/2 one finds, in the heavy quark limit 

/ D ^/ D »^(474±14)MeV (55) 
/ Dl/2 ^(417±13)MeV. (56) 

These values for fo, Jd* are much larger than the estimations given by lattice QCD 
even adding a 10 % error due to quenching (used in Appendix B) : 

f D = (216 ± 36) MeV 

f D * = (258 ± 52) MeV . (57) 
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For /d(o+) one finds, in lattice QCD, keeping only the statistical error [3"4] . 

/ D(0+) = (122 ±43) MeV . (58) 

This latter value is consistent with the value obtained in the quark model of Veseli 
and Dunietz ^H] : 

/ D(0+) = (139 ±30) MeV . (59) 

These theoretical estimations of the /d 1/2 or fo(o+) decay constants, that become 
equal in the heavy quark limit, are not homogeneous in their methods. To make 
the panorama somewhat clearer, we summarize the results in Tables 2 and 3. In 
Table 2 we give the results of the different methods at finite mass, together with 
the phenomenological determination of the present paper. In Table 3 we give the 
results of the methods in the heavy quark limit, dividing the invariant y/mqj 1/2 by 
y /m D{o+) with m D ( +) = 2290 MeV of the Belle experiment (Appendix A). Although, 
of course, this choice is somewhat arbitrary, one can thus qualitatively compare with 
the finite mass results. 



Theoretical method 


/d(o+) 


QCD Sum Rules 


(170 ±20) MeV 


Lattice QCD [H] 


(122 ±43) MeV 


Veseli-Dunietz quark model [TH] 


(139 ±30) MeV 


Present phenomenological 
determination from B — > D**tt 


(206 ± 120) MeV 



Table 2. Theoretical predictions for the decay constant /d(o+) i n the different 
methods, at finite mass, compared with the phenomenological determination of the 
present paper. 
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Theoretical method 


V • D (° + ) 


QCD Sum Rules [25] 


(304 ± 40) MeV 


QCD Sum Rules j3UJ 


(377 ± 53) MeV 


QCD Sum Rules 


(238 ± 66) MeV 


(correcting for Bit continuum) [3Ti] 




Bakamjian-Thomas quark model |T7j 


(417 ± 13) MeV 



Table 3. Theoretical predictions for y/rriQ fy 2 in the heavy quark limit, divided by 
y/ m D(o+), with m D ( +) = 2290 MeV of the Belle experiment. 

The table shows a very scattered set of results, but there is the general trend 
that the decay constant is much larger in the methods that use the heavy quark 
limit. The largest value is obtained by the BT models. The subleading correction is 
negative. Also, we observe that the phenomenological determination of the present 
paper, that has a large error, agrees within errors with the methods including finite 
mass corrections. 



4 Predictions for Class II decays. 

Let us finish our discussion giving predictions for the rates of the color suppressed 
decays, using the range (|4*5Jl for the decay constant /d 1/2 - From the rates 



"<2 „ ^2 



2 



n l/2 n\ _ Gp , |2 P _P 1 I" _ m <• f*B(2 



2 



v y 47r m D 2 1 ' J 

we obtain the branching ratios 

BR (B° -> Do /2 V) = (2.8 ± 2.0) x 10" 4 



BR (F -> Z^ 2 V) = (2.2 ± 1.5) x 10- 4 . (61) 



The central values are large enough that could in principle be measured. These rates 
are independent of the IW function Ti/ 2 (u>), while they depend on the non- vanishing 
decay constant fz> 1/2 . 
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Heavy quark symmetry plus factorization predicts 

BR (B° -> B-T V) = BR (b° _> £>f 2 °7r°) = (62) 

because of the vanishing of the decay constants (|27)1. However, it is worth 
noticing that, because of the large experimental errors and theoretical uncertainties 
(spin counting, etc.), from the BR (|2l)j) and using the isospin relation (0) we can 
only have a rather loose upper bound 

B R (B° -> DT V) < 4 x 10- 4 

BR (B° -> Df 2 7t°) < 3 x 1(T 4 . (63) 



5 The rate to excited states in semileptonic B 
decays. 



The values of the functions 71/2(10), 73/2(10) at w = 1 and their w-dependence 
gives predictions for the semileptonic (SL) decay B — > D**£u branching ratios in the 
heavy quark limit. The differential decay rates for B — > £)j (j = |, |)&/ write [T7] 



dw 



^f\V cb \ 2 2r%w + l)(w 2 -l)^ 
w + l)(l- r) 2 + 3iu(l + r 2 - 2rw) 



T3/2I 



10 



dw 



^f\V cb \ 2 2r\w + l){w 2 -lf/ 2 
(w - 1)(1 + r) 2 + w(l + r 2 - 2rw) 



rfr (b -> Di 72 ^// 

dw 



dV (B -> Dq /2 £z/ 
dw 



^\V cb \ 2 Ar\w-l)(w 2 -lf 2 



(w - 1)(1 + r) 2 + 4iu(l + r 2 - 2rw) r 1/2 (w) (64) 



^|^| 2 4rV 2 -l) 3/2 (l-r) 2 |r 1/2 H 



where r = ^ and D denotes the corresponding Z)j meson. 
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For completeness, we write down the corrresponding formulas for the ground 
state : 

dT{B -> Div) = G|r4 |K&|V(w2 _ 1)3 / 2(1 +r)2|eH|2 



dw 

dT(B -> D*tv 
dw 



48tH 



^f|K b |V(i + -)(^ 2 -i) 1/2 

407T 5 



(w + l)(l-r) 2 + 4w(l + r 2 -2rw) |^(w)| 2 . (65) 



The situation is given in Table 4, a slight modification of the predictions of ref. [T7] . 



Semileptonic mode 


Experiment 


Model 




(2.14 ±0.20) x 10~ 2 


(1.95 ±0.45) x 10" 2 


B -> D*£v 


(5.44 ±0.23) x 10~ 2 


(5.90 ± 1.10) x 10" 2 


B -> B 3 / 2 ^ 


(a) (2.4 ± 1.1) x 10~ 3 

(b) (4.4 ±2.4) x 10~ 3 

(c) (3.0 ±3.4) x 10~ 3 


(6.3t 3 ;°) x io- 3 


B -> L>i /2 £z/ 


(a) (7.0 ± 1.6) x 10- 3 

(b) (6.7 ±2.1) x 10~ 3 

(c) (5.6 ± 1.6) x 10~ 3 


(4.o±i;l) x io~ 3 


B -> Z>} /2 &/ 


(b) (2.3 ±0.7) x 10~ 2 
"wide" D** -»• (D + D*)tt 


(6 ±2) x 10" 4 


B -> Bq /2 £z/ 


(6 ±2) x 10" 4 



Table 4. Comparison between rates for B — > Dj£v decays and the model described 
in the text. The data are from (a) ALEPH |3S| |3S|, (b) DELPHI |H| [37| [HE] and 
(c) CLEO [H] experiments. For the elastic IW function we adopt the values (B.5), 
(B.6) of Appendix B. "Wide" stands for unidentified (D ± D*)ir events forming a 
wide bump. In the text we discuss a new DELPHI analysis [35] and very recent 
data from the Belle collaboration 140 . 



To make predictions for the SL rates we need an input on the IW functions 
T i/2( w )> t 3 / 2 (w). First, we must take into account that it is reasonable to expect 
that the n = IW functions give a sizeable contribution to Bjorken and Uraltsev 
SR. Making this assumption, Uraltsev SR is very constraining on the difference 



|r 3/2 (l)P 



l T l/2 



1)| 2 , that should be not far away from -r. This is the case for the 



BT model values ()21|) . and, more importantly, consistent with the values that we 
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have found from non-leptonic decays (J2*U|) using the Belle data. As an example, 
we will then adopt the values |r 3 / 2 (l)| and |ri/ 2 (l)| given by the BT model (J2T|) 
and allow nevertheless a ±50% departure for the values of the slopes. The lower 
values of the slopes would be in agreement with ref. [27] . This gives the range 
of model predictions for SL rates in the table. Although the B — > D J £v rates 

1/2 

are in reasonable agreement, for the B — > DJ £v rates there is a problem. We 
find BR[B -> (£>o /2 + I>} /2 )&>] = (1.2 ± 0.4) x 10" 3 . The DELPHI experiment 
gives, once the BR[B — > {D^ + D^ 2 )£u] is subtracted, a large branching ratio of B 
decays into "wide" D** mesons decaying into (D + D*)tt. We will call this branching 
ratio BR^ PHI \B -> (D + L>*)tt£i/] = (2.3 ± 0.7) x 10~ 2 , that we report in the 
last line of Table 4. On the one hand, this BR is one order of magnitude larger 
than our prediction for BR[B — > (D^ 2 + D\^ 2 )£v\. On the other hand, keeping 
only to the DELPHI experiment, one obtains the sum BR[B — > (D + D*)iu] + 
BR[B -> (D 3 / 2 + D 3/2 )£z/] + BR wide [B —> (D + D*)tt£v] = (11 ± 1.6) x 10' 2 , 
that already saturates within errors, that are however large, the total semileptonic 
width BR[B ^ £u + anything] = (10.73 ± 0.28) x 10~ 2 . There are therefore two 

1 /2 1 /2 

problems. On the one hand, other wide states besides D ' +-D 1 / have to contribute 
to BR wi( ie[B — » (D + D*)7t£u}. These could be radial excitations or higher orbital 
excitations, that are in principle allowed due to the large phase space available. 
The experimental width is so large that it includes high masses. Moreover, the 
multiplicity of higher excitations grows with the mass. On the other hand, it is 
curious that considering only the modes B — > (D + D*)£i>, B — > (D + D*)it£v the 
total semileptonic width is already saturated, and one could wonder why there is 
no place for decays into multipion modes B — » (D + D*) + nrr (n > 1) and why 
they are not observed. Presumably, due to phase space, these could not come from 
modes of the type Dp but could come from various D**n. 

A recent new analysis by DELPHI [SH] confirms and makes more precise this 
situation. The total branching ratio into D** (narrow and broad) is measured to be 



with the decay final states dominated by the D(D*)ir channels. The dominant 
contributing channel is a broad state decaying into D*ir, i.e. a state D** or J p = 1 + 
with a mass M = 2445 ± 34 ± 10 MeV and a width T = 234 ± 74 ± 25 MeV. On the 




(66) 
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other hand, broad states favor a production with a maximum close to threshold. 
On the other hand, DELPHI bounds the branching ratios into Dim and D*7ttt final 
states. 

On the other hand, very interesting recent data have been published by Belle at 
the last Lepton-Photon Conference |40| . that gives the following branching ratios 



BR (B- -> D + n~rv) = (0.54 ± 0.07 ± 0.07 ± 0.06) x 10 
BR (B- -> D* + TT-rv) = (0.67 ± 0.11 ± 0.09 ± 0.03) x 10 
BR (b° -> D G n + rv) = (0.33 ± 0.06 ± 0.06 ± 0.03) x 10~ 2 



2 

-2 



BR (5° -> D*°7r + ri7) = (0.65 ± 0.12 ± 0.08 ± 0.05) x 10~ 2 . (67) 

These values are to be compared with the theoretical expectations of Table 4. We 
find, considering only the central values of the model and taking into account the 
relevant branching fractions of the different Dj computed in Section 2, 

BR(b~ -> D + ir~£~v) = BR (b° -> D°7r + rv) = 0.43 x 10~ 2 



BR (B- -> D* + v-rv) = BR(b° -> D*°7r + ri7) = 0.34 x 10~ 2 (68 



In view of the uncertainties, there is a fair agreement between our predictions 
and the Belle data. 

In conclusion, there seems to be a potential problem concerning the DELPHI 
semi-leptonic data on B —> (D^ 2 +D^ 2 )£i' decays that should be addressed in future 
experiments. If BR wic i e [B — > (D + D*)^^} had to be attributed to D\' 2 -\-D\ , then 
71/2(1) would be much larger than T3/2(l), in contradiction with the expectations of 
Uraltsev SR. This is at odds with the Belle non-leptonic data studied in the present 
paper. The study of the Dj wide states is not an easy experimental task. A recent 
Tevatron DO experiment sees clearly in SL B decays the narrow states j — | but 
has not given a measurement of the wide ones j ' — ~ [32] • On the other hand, we 
find agreement with the very recent Belle data on B — > D(D*)7r£v. 

6 Conclusion. 

In conclusion, we have shown within a simple factorization model, tested in 
the well-measured B — > D(D*)tt decays, that one can extract information on the 
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Isgur-Wise functions at zero recoil 71/2(1), 73/2(1) from non-leptonic data on B — ► 
D** + 7i~ (Class I). Combining with B~ — > D**°tt~ (Class III) one can obtain the 
non- vanishing decay constant /d 1/2 of D**(j = \). Special care has been taken in 
the determination of the interference sign between the 7r and D** emission diagrams. 
The ranges obtained for 7i/ 2 (l), t 3 / 2 (1) are consistent for both types of modes, with 
the expectations of Bjorken and Uraltsev sum rules, and with the predictions of 
the Bakamjian-Thomas quark model of form factors. Moreover, the range of values 
found for the decay constant /d 1/2 agree with most theoretical expectations. We 
predict sizeable rates of the Class II decays B° — > D**°7r° that could be measured 
in the near future. On the contrary, for D**(j = |), Class II decays should be 
suppressed due to the vanishing of the decay constant f 3 / 2 in the heavy quark limit. 

We must warn that 1/wq corrections could be large and upset the results of 
the present stage of this analysis, as discussed in Appendix D. Also, we point out 
a problem with present DELPHI data of semileptonic decays for the total rate to 
excited states. Very recent Belle data on B — ► D^iriv seem in good agreement 
with the theoretical expectations. 
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Appendix A. Belle data on B —> D**7r decays. 

For the sake of clarity on how we extract the branching ratios of the different 
B — > D**7i decay modes, and how we handle the errors in the text, we reproduce 
here the Belle data on Class I [3] and Class III decays jlj. 
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M 2 3/2 




(2459.5 ± 2.3 ± 0.7l£§ 


) MeV 




I 


,3/2 
2 






(48.9 ± 5.4 ± 4.2 ± 1.0) MeV 


B 


(B° - + 7T- 


)* 






(3.08 ± 0.33 ± 0.09lg;ra 


) x 10" 4 


B 








) 


(2.45 ± 0.42l°;| ±g^?) 


x 10" 4 


Ml' 2 


(2428.2 ± 2.9 ± 1.6 ± 0.6) MeV 




r 


,3/2 
1 






(34.9 ± 6.6t 4 ;£ ±4.l) 


MeV 


B 


- dT v-) 


B[ 


+ ^ 


) 


( 


>.68±0.6018:S + -° g) 


x 10~ 4 


Ml 12 


(2290 ± 22 ± 20) MeV 


r l/2 
1 


(26.7 ±3.1 ±2.2) MeV 


B 




)B 






< 1.2 x 10" 4 


Ml 12 


(2428 ± 2.9 ± 1.6 ± 0.6) MeV 


r l/2 


(380 ± 100 ± 100) MeV 


B 




B 




) 


< 0.7 x 10" 4 



Table 5. Data of the Belle Collaboration for the masses, widths and branching 
ratios to D** states D j 7 for Class I decays B° — > D** + tt~ jS]. 



M 2 3/2 


(2461.6 ± 2.1 ± 0.5 ± 3.3) MeV 


F 3/2 
1 2 


(45.6 ± 4.4 ± 6.5 ± 1.6) MeV 


B 


(b- -> dT V- 


)* 






(3.4 ±0.3 ±0.6 ±0.4) x 10" 4 


B( 


'b- - d 3/2 V) 




^ 


) 


(1.8 ±0.3 ±0.3 ±0.2) x 10" 4 


M 3/2 


(2421.4 ± 1.5 ± 0.4 ± 0.8) MeV 


rf 2 


(23.7 ±2.7 ±0.2 ±4.0) MeV 


B{ 


'B- _> uj/a V) 




^3/2 ^ 


) 


(6.8 ±0.7 ±1.3 ±0.3) x 10" 4 


Ml 12 


(2308 ± 17 ± 15 ± 28) MeV 


r l/2 
1 


(276 ± 21 ± 18 ± 60) MeV 


B 


(5- - z^ /2 v- 


)B 


(^ /2 ° - D + 7T-) 




(6.1 ±0.6 ±0.9 ± 1.6) x 10" 4 


Ml' 2 


(2427.0 ± 26 ± 20 ± 15) MeV 


if 2 


(384±^ 7 ± 24 ± 70) MeV 


B{ 


- Di /2 v-) 


B{ 


^V 2 o _> D * +7r ^ 


) 


(5.0 ±0.4 ±1.0 ±0.4) x 10" 4 



Table 6. Data of the Belle Collaboration for the masses, widths and branching 
ratios to D** states D j j for Class III decays B~ -> D** 7r~ 
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Appendix B. Testing the factorization model in 
B D(D*)tt. 

In this Appendix, in order to check qualitatively our simple factorization model 
applied to B — > D**tt decay, we use it to describe the well measured decays into the 
ground state B — > D(D*)tt. The data of the PDG ^H] is given in Table 6, together 
with our predictions that follow from the following simple formulae. 
From the definitions 

< n-(q)\A^\0 > = f v q fl 
< D°{p')\A ll \B-ip) > = ^/m D m B £(w) (v + u% 



v)v, 



<D°(p')\A li \0>=f D p' fi 
<D*°(p',e)\V,\0> = f D m D , ^ 

< n-( q )\V,\B-(p) > = (p» + q»- T ^^v')j ff( P ' 2 ) + 

P P (b.i; 

we obtain the matrix elements, satisfying the isospin relation (0), 

-\ r— — — — „, \ /„,, i iW £/„.. \ i „ _2 t f-KBf2 \ 



A[B -> D 7T J = a ly /m D m B (m B - m D ) (w + l)f^(w) + a 2 m B f D f£ (m D ) 
A (b° -> D+tt") = a ly /m D m B (m B - m D ) (w + l)f n £(w) 

A(b°^ Z>V) =-a 2 -Lmlf D fZ B (ml) (B.2) 







)= P [ 









a lv /m D *m B (m B + m D *) f^(w) 
+ a 2 2m D *m B f D . f+ B {m 2 D * 

A (b° -> L>* + 7r~) = ax — — ^m D *m B (m B + m D *) f^(w) 
v 7 mo* 



A l B " _> D*°7r° = - a 2 -= 2m D ,m B /™K>*) • (B.3) 
v 7 V2 m D * 

For the decay constants /©, fp* we use the values of lattice calculations within 
the quenched aproximation f D = 216(H)(5) MeV |32j, f D , = 258 (14) (6) MeV 
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|33| . where the first error is statistical and the second is systematic. Assuming a 
10 % uncertainty due to the quenching approximation and adding all the errors in 
quadrature, we adopt the values 

f D = (0.216 ±0.036) GeV 

fjj* = (0.258 ± 0.052) GeV (B.4) 

For the form factors /o B (g 2 ), f+ B {q 2 ) we use (B2Jl - (|33jl . an d f° r the Isgur-Wise 
function, we use the parametrization given by the BT model (last reference of [Tl]). 




(B.5) 



that we have used in [13] to fit Belle data on B -> D*£v |Sj, that gives F*(l)\V cb \ = 
0.036 ± 0.002 and p 2 = 1.15 ± 0.18. In conclusion, from .P(l) = 0.91, we adopt the 
ranges 

\V cb \ = 0.040 ±0.002 

p 2 = 1.15 ±0.18 . (B.6) 

We add the theoretical errors in quadrature, that gives, in amplitude, a 10 % error 
for Class I decays and a 30 % error for Class II decays. Adopting the values 

ai = l , a 2 = 0.2 (B.7) 

obtained by the perturbative calculations [22]; and considering the uncertainties 
given in (B.4) and (B.6), we obtain the predictions of Table 7. 
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Modes 


Experiment 


Factorization 


B D(D*)tt 




a x = 1, a 2 = 0.2 


BR(b° -> 


(2.76 ± 0.25) x 10~ 3 


(3.4 ± 0.7) x 10~ 3 


BR(B ^ D°7T ) 


(2.7 ± 0.8) x 10 -4 


(0.6 ± 0.4) x 10 -4 


BR(B~ -> D°7T-) 


(4.98 ± 0.29) x 10" 3 


(5.2 ± 1.0) x 10" 3 




(2.76 ±0.21) x 10~ 3 


(3.3 ±0.7) x 10" 3 




(2.7 ±0.5) x 10" 4 


(0.8 ±0.5) x 10" 4 




(4.6 ±0.4) x 10" 3 


(5.3 ± 1.0) x 10" 3 



Table 7. Data on the branching ratios of B -> D(D*)tt from PDG 2004 [EJ, 
and the predictions of the factorization model with the perturbative values a\ = 1, 
a 2 = 0.2. The errors in the predictions come from the uncertainty on the decay 
constants of D(D*) mesons, the value of \V c b\, the Isgur-Wise function and the form 
factors fl B (q 2 ), fo B {q 2 ) AH"®. The theoretical errors are added in quadrature. 

A first remark on the experimental data of Table 7 is that the isospin relation (0) is 
roughly satisfied by real amplitudes, i.e. without the need of FSI phases, like in the 
naive factorization model. However, it is clear from this table that Class II decays 
are underestimated using the perturbative value a 2 — 0.2. There is no theoretical 
reason, unlike the case of Class I decays, in which there is emission of the light meson 
7r [7], to have approximate factorization in the case of D(D*) emission. There are 
non-perturbative corrections to factorization, as pointed out in Appendix C, that 
suggest an effective value for a 2 - In Table 8 we give the results for 

ai ^ 1 , a 2 = 0.3 (B.8) 

that agree with the data within errors. 
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Modes 
B D(D*)tt 


Experiment 


Factorization 
a x = 1, a 2 = 0.3 


BR(b° -> 


(2.76 ± 0.25) x 10~ 3 


(3.4 ± 0.7) x 10~ 3 


BR(B ^ D°7T ) 


(2.7 ± 0.8) x 10 -4 


(1.3 ± 0.9) x 10 -4 


BR(B~ -> D°7T-) 


(4.98 ± 0.29) x 10" 3 


(6.0 ± 1.2) x 10~ 3 


55 (B° -> D*+7T-) 


(2.76 ±0.21) x 10~ 3 


(3.3 ±0.7) x 10" 3 


55 (5° -> 5*°7r°) 


(2.7 ±0.5) x 10" 4 


(1.7 ± 1.1) x 10" 4 


5/2(5- -> 5*°7T-) 


(4.6 ±0.4) x 10" 3 


(6.3 ± 1.3) x 10" 3 



Table 8. Same as Table 6 with the effective values a\ = 1, a 2 =■ 0.3. 



Now Class II decays are in better agreement and the overall picture seems reasonable. 
In the estimation of the 5 — > D**tt decays in Section 3 we adopt the values (B.8) 
for a\ and a 2 . 



Appendix C. Corrections to factorization in 
B — ► D**7T decays. 



In the simple-minded factorization approach one adopts the perturbative Wilson 
coefficients a± = 1, a 2 — 0.2. Moreover, in this approach, the amplitudes of Class II 
decays to j — | states vanish in the heavy quark limit A(B° — > D^n ) = due to 
the vanishing of the decay constant / 3 / 2 - 

Let us here discuss how the analysis would be modified by taking into account 
non-perturbative corrections to factorization, following Neubert [3]. In ref. [3] only 
the decays 5 — > D^* 2 tt are discussed. We extend the formalism to 5 — > Dy 2 7i. 

In the case of the j = \ states D** emission is allowed. Therefore, we write the 
amplitudes of Class I and Class II decays 5 — > D**, 2 tt ( J = 0, 1) in the form 



A (5° 



A (5 



D 



D 



1/2 



71 



a 



.4 



fact 



-ftO 



5 



1/2 O n 
J 71 



a/ M2 A fact (5° 



n l/2 + - 
Dj 71 



n l/2 + - 
DJ 71 



(c.i; 



with 



e//,l/2 



Cl{fj) + 



l /2 (n)] ±c 2 (/i) eJ /2 (//) 
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„e//,l/2 



c 2 (/x) + 



1 + £i /2 (/x) +ci(/x) 



1/2, 



(C.2) 



where the hadronic parameters e 1 / 2 , e^ 2 describing the non-factorizable contribu- 



tions are given by the matrix elements 

< D 1 / 2 + 7i-\(du)(cb)\B° > = \l + e\ /2 (i2) 



A 



fact 



D)' 2 V 



< D 1 ] 2 + n-\(dt a u)(ct a b)\B° > = \ el /2 (fi) A fact (b° -> D 1 ] 2 



71 



(C.3) 



We make explicit the upper script j = |, since the situation is quite different for 
the decays into j = | states. Following pQ , we consider the large- N c counting rules 



d = 1 + O (l/iV c 2 ) 
ei = O (l/iV c 2 ) 



c 2 = O (l/iV c ) 

e B = 0(l/N c ) . 
Keeping the terms up to order l/N c included, one finds 

a\ f /,1/2 (/i) = of r *(/i) = 1 

a^ 1/2 (/i) = af rt (/i) + ^s /2 (/i) = 0.2 + 4 /2 (/,) 



(C.4) 



(C.5) 



The departures relative to the naive approximation presented above are given by the 
non-perturbative coefficient s^ 2 (fi). These quantities do not affect Class I decays, 
but only Class II and Class III [cf. the isospin relation (1)]. 

In the analysis of Appendix B on the well measured B — > D{D*)n decays, we 
have found a\ = 1 and = 0.3. Although there is no firm theoretical argument, 
we have adopted in the text the same values, i.e., 



a 



e//,l/2 



e//,l/2 
u 2 



0.3 . 



(C.6) 



Going now to the case of j — § states, we define 



A (B° -» Dj 



3/2 



-V2 (B° - DT u 7r u j = sT A fact (£" - 
where, due to (C.4), 



7T 



,3/2 On 



a 



e//,3/2 



-ftO 



.3/2 



,3/2 



(C.7) 



_e//,3/2 



Ci(/i) + 



C2M 



1 + ef (/i) + c 2 (/i)4 /Z (/i) = «f rt (/i) = 1 • (C 



.3/2 



pert 1 
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The hadronic coefficients e^ 2 (fi), e^ 2 (fi) are defined like in (C.3). We have kept 
the notation e% 2 in (C.7) because, switching off non-perturbative corrections, this 
coefficient does not have as a limit a non- vanishing perturbative coefficient, unlike 



a 



e//,3/2 



. Since the amplitude A(B° — > D 



7T 



vanishes in the heavy quark limit, 



because / 3 / 2 = 0, the amplitude chosen for the normalization is the one that is 
allowed, the ir~ emission one. 

Taking into account a non- vanishing coefficient e^ 2 in (C.7), 



r(# 



dT V) = r (5° - £>i 



3/2 0^0 



16tt 



= 1 
2 



.3/2 



[1 + rY 



1 6r 3 



T"3/2 



1 + 

2r 



(C.9) 



where r = 

' l+r 2 



3/2 



mi" = 0.46. From the upper limits (fMJ) and the central value (fTH|) for 



T"3/2 



2r 



we can infer an upper limit for |£g^ 2 |, namely 

< 0.90 . 



I 3/2, 
I 



(CIO) 



As expected, since the upper bounds (}63|) are rather loose, we obtain a large upper 

3/2 i 

bound on \e R . 



Appendix D. Remarks on 1/niQ corrections. 

Using the formalism of [28: and assuming factorization, one can in principle 
compute the analytical expressions of the I/tuq (Q = b or c) corrections to the 
rates B — > D^ix. Let us consider Class I decays, 5° — > D^ix~ . Many subleading 
form factors contribute and, although a theoretical effort has been made in their 
estimation for the j ' — | states within the QCD Sum Rules approach [45) . we do 
not have presently at our disposal an estimation of all the subleading form factors 
defined in j2Bj- Therefore, we are not able at present to make an estimation of these 
corrections. However, a formal expansion can be done for these decays to pions, and 
subleading quantities can be estimated in some approximation, as we explain now. 

Let us consider the most important contributions at w — w max = Wq, that 
corresponds to q 2 = 0, the value for pion decays : 

w = w = £ = —— = 1.3 (D.4) 

IrriBmr) 2r 
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where is the mass of the corresponding Dj meson and r = mjj/iriB- Therefore, 
one can express all the mass factors in terms of wq and a common overall scale. 

Then, for ir decays, using (D.4), there are two small parameters that characterize 
the corrections to the rates, namely 

w - 1 = 0.3 and (Q = c,b) . (D.7) 

It is therefore convenient to classify the subleading corrections to the rates as being 
of successive orders 

( »o-ir(^)'. to,. 

We decide to retain only the subleading orders contributing to the rate (t = 1) with 
the dominant order in (u? — 1)> namely s = —1/2. 

Using the formulas of ref. j2Ej, neglecting higher orders of the type (D.8), and 
keeeping the dominant order (wq — l)" 1 / 2 ( ) ? we fi n d 

r (B° - D?*-) = T (B° - (l + 55) (D.9) 



with 



bT = o 

2v/2 A£ 3 /2 1 



.3/2 _ ^V^ ^^3/2 
01 



— 1 2m c 
i V2A£ 1/2 ( 3 1 



v^o — 1 \2rrib 2m c 



s ! /2 _ 3V2 AEi/2 / 1 + 1 \ (D 10) 

" — 1 \2m b 2m c J 

and T denotes the leading rate. Numerically, some of these terms of order I/tuq 
are not small, as they are respectively of the order 

S 3 2 /2 = ^ /2 ^0.7 5] /2 ^0 5d /2 ^1.3 (D.ll) 

for AE 3/2 ^ AE 1/2 = 0.4 GeV, m c 1.5 GeV, m b ^ 4.8 GeV and w - 1 = 0.3. 

For j ' — I states, the trend is not in the right direction to explain the different 
central values (fTTj) . This gives an idea of the type of uncertainties induced by the 
1/rriQ corrections, that are large. One can guess that the extraction of |ti/ 2 (wo)| 2 , 
\t 3 / 2 (wq) | 2 made in Section 2 is uncertain by about a factor 2 due to these corrections. 
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Therefore one could assume a reasonable additional 40% uncertainty on the values 
(HHJ) given for |ti/ 2 (wo)|, 1 73/2 Oo) I- 

Although these estimations give large corrections, this impression could be wrong, 
since many subleading form factors contribute |28j . and we do not know the mag- 
nitude or sign of the neglected terms. Moreover, we have taken only the leading 
order in the expansion (D.8). The aim of this exercise has been to emphasize that 
the corrections in 1/rriQ are possibly large. If this was actually the case, this would 
upset the results of the present stage of this analysis. 
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